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sodium hydrate. Acids were also used for dissolving the metallic splash, 
with negative results. 

If the xenon is present as a metallic compound, it must therefore go into 
solution as a compound on dissolution. Silicon hydrides, when treated with 
sodium hydrate, gives sodium silicate ; it may be that xenon behaves in a 
similar manner, and that a xenate of sodium is formed. 



The Supply of Energy from and to Atmospheric Eddies. 

By Lewis F. Eichaedson. 

(Communicated by Sir Napier Shaw, F.R.S. Received March 9, 1920.) 

Introduction, 

This paper extends Osborne Reynold's theory of the Criterion of 
Turbulence, to make it apply to the case in which work is done by the 
eddies, acting as thermodynamic engines in a gravitating atmosphere. For 
simplicity, the air is supposed to be dry. The gist of the solution is con- 
tained in equations (4*16) and (5*3) below. 

Some analogies to thermodynamics are then introduced, including the con- 
ception of an eddy-heat-per-mass, and its diminution when the air expands. 
These conceptions are helpful in discussing the manner in which eddies 
diffuse. I am indebted to Mr. W. H. Dines, F.R.S., for his valuable advice 
at many stages of this work. 

While it was being finished off, a treatment of the Criterion of Turbulence, 
by Wilhelm Schmidt, of Vienna,* has come to hand. For the supply of 
energy from the wind, he gives much fuller observational data than I have 
done. But, for the thermodynamic loss, he concludes that the activity per 
volume has the average value hdW/d¥ where h is the gas constant in 
erg units, 6 is the absolute temperature, and h the height. This is incom- 
patible with the result arrived at in equation (4*15) below, and is, moreover, 
difiioult to reconcile with observation. I am* inclined to think that 
Prof. Schmidt has accidentally admitted a factor, p^"^^, in the diffusion 
equation. On integrating partially, and taking account of boundary con- 
ditions, this change would explain the discrepancy, if the turbulence were 
independent of height, as he assumes. 

^ * Ann. der Hydrog. ii. M^r. Meteorol,,' November- December, 1918. 
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{Added May 25, 1920. — At the meeting at which this paper was discussed 
Mr. 6. I. Taylor mentioned some studies of stability which formed part of 
his Adams Prize Essay in 1914. The oscillations of superposed laminae of 
different density became unstable when a criterion was reached which 
resembled (5*4) below, except that the numerical coefficient contained a factor 
of 2 or of 4 depending on the number of the laminae. By a rough energy- 
method he had arrived at a numerical coefficient agreeing with that of (5*4). 
It is to be hoped that these investigations will appear in print.] 

I. The Available Forms of JEnergy. 

Since energy is conserved, we might find the supply to eddies by counting 
up the losses from all the other forms. 

In order to relieve ourselves of the trouble of considering effects at 
vertical walls, let us consider a large volume of atmosphere, say that 
standing on a land area measuring 10 km. x 10 km. or larger. Let this area 
be A cm.^. 

Let p^ be the mean pressure at the earth's surface. Then the mass 
standing upon A is Ap(^/g = M, say. Let dM. be an element of this mass. 
The changes of the energy present in this mass may be classified as 

(i) Change of. intrinsic thermal energy AI, 

AI = |7^ . A(9 . dM, (1-1) 

where 7^ is the thermal capacity per mass at constant volume, and is fhe 
temperature. 

(ii) Change of gravitational energy AT, 

Ar=:A|.^MM, (1-2) 

wheie h is the height of dM. and g is the acceleration of gravity. 

(iii) Imagine the wind velocity to be smoothed, so as to remove the gusts. 
Let Vx, %, ^H be the actual winds ; let Vx, v^, v^ be the means of the same. 
The means may be supposed to be taken over a time-interval, long compared 
with gusts, but short compared with the passage of barometric depressions. 

Denote deviations from the mean by dashes, so that 

'2^x = % + < ; "0^ = ^Y + '^y' ; % = '^H + V, (1-3) 

Square, add and multiply by half the density p, and again take the mean, 
denoted by a bar, over a similar interval of time. Then, as is well known, 
the terms containing products of means and deviations vanish, leaving 

so that the kinetic energy may be divided into that associated with the 
mean and that associated with the deviations. Denote these parts, for our 
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large atmospheric block, by E and E' respectively, and their increases by 
dE and AE'. The object of this paper is to study AE'. 
By the continuity of energy 

(AI + Ar+AEH-AEO/A^ = G + E+W, (1-5) 

where G is the rate at which heat-energy is flowing into the large portion 
of atmosphere from the vegetation and ground or from the sea. And E 
is the excess of the radiation flowing in over that flowing out. 

And W is the rate at which work is being transmitted across the vertical 
walls in the form of (pressure) x (velocity). As the size of the area A 
increases, W dwindles relatively to the other terms. Let us neglect W. 

The immediate effect of the radiation E is to alter the temperature of the 
air, vegetation, land or sea. Afterwards, these temperature-changes may 
cause eddies. But in this paper we suppose that the temperature distribu- 
tion is known. Therefore we may ignore radiation, provided that we also 
ignore the changes of temperature which it is at the instant producing. 
Let Aol be the change in the intrinsic energy not due to radiation. Accord- 
ingly (I'o), becomes 

(Aol + Ar + AE + AE')/ A^ = G. (1-6) 

This equation may be further simplified, for Mr. W. H. Dines has shown 
that there is a connection between the changes of gravitational and intrinsic 
energy.* On inserting the appropriate constants, 7^ the thermal capacity 
per mass at constant volume, and h the " gas constant '' defined by the 
characteristic equation 

23 = hpO^ (1**7) 

where p is density, Mr. Dines' theorem takes the following form: ''In a 
column of clear air hounded hy vertical walls, and extending from the ground 
upwards to a height at which the pressure is negligible, the change in total 
gravitational energy, reckoned froon the ground, is hj^^ times the corresponding 
change in the total intrinsic energy!' (1'8) 

Now, ^1, + 6 = 7^, where 7^ is the thermal capacity per mass at constant 
pressure. The theorem (1*8) may therefore be interpreted thus : a column of 
air, of unit mass, in cooling, gives out energy as if it had a specific heat jp. 
Of this energy, a part, 6/7^, is due to the sinking of the column. (1'9) 

So equation (1*6) simplifies to 

l^^+^+^r^G. (1-10) 

7t, at at at 

The whole energy is sometimes correctly classified in ways involving the 

* ' Q. J. lioy. Met. Soc.,' p. 188, July, 1913, 
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pressure p, Margules,* for instance, calculates the energy associated with 
local irregularities of pressure. In the present scheme this energy comes in 
(i) and (ii), and, to include pressural energy in addition, would be to 
count some of the energy twice over. It seems to be most natural to 
regard pressure only as concerned in the transfer across surfaces. Once 
across, energy appears in some other form — gravitational, intrinsic, or molar 
kinetic. 

II. The Siif^ly from the Mean Wind, 

Osborne Eeynoldsf showed that the kinetic energy of eddying was 
increased by the work done by the eddy-stresses upon the corresponding 
rates of mean strain. 

In the atmosphere, high above trees and houses, many of the rates of mean 
strain are negligible, and the activity may be abbreviated to 

activity per volume = xh -—-{•yh -^y , (2'1) 

oh ' vh 

where xh, yh, are the shearing stress on horizontal planes and where h is the 
upward co-ordinate. 

If fji be the " eddy viscosity " defined as 

__ eddy shearing stress <'9-9\ 

corresponding rate of mean shear' 

Then equation (2*1) may be written 

activity per volume = /^xh ( g^ j + a^yh f g^ j , (2-3) 

where />6xh and /Xyh are possibly not equal, as suggested by the author,^ 
expressing the rate at which the mean-wind gives up its energy to eddies. 
ISTumerical data, based on observation, are collected by H. U. Sverdrup.§ 

III. The Loss hy Molecular Viscosity, 

Reynoldsll balances the above gain of energy against the energy converted 
into heat at a rate per volume equal to 

* " The Mechanical Equivalent of Pressure," ' Ahbe's Translations,' 3rd series, 
Smithsonian Publications. 

t Lamb, * Hydrodynamics,' 4th ed., §369. 
X ' Phil. Trans.,' A, p. 1 (1920). 
§ * Geo. Inst., Leipzig,' 2nd ser.. Band II, Heft 4. 
Lamb, 'Hydrodynamics,' 4th ed., p. 369. 
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where fjbc is the molecular viscosity and dashes denote deviations from the 
mean. 

This rate of loss will vanish if the eddying energy vanishes, for then 
%', %^ %' will be zero everywhere. Further, it is clear that the loss by 
molecular viscosity increases with the eddying velocities and with the 
stnallness of the eddies. 

The expression (3'1) when integrated over the volume standing on the 

larffe area A will be denoted by—- (-ts— ) . 

The suffix is necessary because it is not the only part of -^. 

To gain some rough idea of the numerical value of the dissipation, let us 
imagine that the complicated irregular motion is analysed into a series of 
regular motions superposed upon one another according to the method of 
Fourier. In other words, let each component of the instantaneous velocity 
be expanded in sines and cosines of each of the co-ordinates. Let us pick out 
for examination the terms 

v^ ~ B sni — cos — ^ sm — 

''^ ""' ''" ; ./ = ^ (3-2) 

, T> nroo wy irli 
q)rr == B cos ~ cos -^ cos — - 



^x ^y ^'H J 



which signifies a motion divided up into rectangular cells, so that the eddies 
would look rather like papers curled up in the pigeonholes of an office desk. 
The size of the cell is n^ by n^^ centimetres in the plane of the motion, and 
71y at right angles to the motion. 

The mean kinetic energy associated with this distribution is got by 
squaring the trigonometrical terms and integrating over a volume of atmo- 
sphere very large compared with the size of the cell. Taking the mean with 
respect to any one co-ordinate replaces the square of a sine or cosine by J. 
So the three successive operations introduce the factor (^)^ and the mean 
kinetic energy is therefore 

lp{v^ '^-^v^'^~{-<^) = i-pM-f + 1} = IpB^ = E7(the large volume). 

(3-3) 

Next to find the rate of dissipation we have to insert v^ \ % ', i^h ' in the 
dissipation function quoted above. The differentiations introduce factors 
tt/tix, 'rr/ny, tt/^h- The term dv^ /dh~j-dv^/dx conveniently vanishes. 
Squaring and taking the mean over a large volume again introduces the 
factor ^. So that the mean rate of dissipation comes to 

2,,,7r^JJ~^ + ^^ + J-\^B-^^^^ ^ (3-4) 

L'^H % ^H J oi (the large volume) 
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The quantity ('^Ix^^ + ^^^+Wh""^) is a measure of the smallness of the 
eddy. I.et us denote it by l~^, so that we may speak of I as the linear size 
of an eddy. If the cell which encloses the eddy were cubical, then I would 
be 0*577 of the edge of the cube. If the cube were drawn out into a very 
long square prism, then I would be 0*707 of the side of the square. If the 
prism were stretched into a plate, I would be the thickness of the plate. 

Then from (3*3) and (3*4) 

|:_2^,,, (3-6) 

Now fie = 1-7 X 10-^ for air at 273° A. 
So if I be say 10 metres 

1 BE' 



W Bt 



:=^ 2*6 X 10-^ see-i. 



So that after 24 hours the eddying energy would have fallen to 0*8 of its 
initial value. The changes which we observe near the ground are much 
more rapid, so that either the eddies must be smaller than the size repre- 
sented by / == 10 metres, or else other causes than molecular viscosity must 
cooperate to destroy the eddying energy. 

We will not now stop to inquire what would be the interaction on the 
average between eddies of different sizes and positions represented by the 
different terms of the Fourier expansion. The equation (3*5) is to be 
understood as a type, not as a complete formulation. 

IV. Convection. 

The rising of cumuli during the daytime and the diminution of turbulence 
in the surface wind on clear nights both shows us that convection has an 
important influence. 

The eddies are themselves the engines by which the available heat is 
converted into their visible motion, so that the rate of conversion must depend 
somehow upon the rate of eddying. Let us explore this question. 

But the eddies are imperfect engines, for they lead to mixing of portions 
of air at different temperatures, and mixing is an irreversible process. If, 
however, we can by any means find the diminution of the intrinsic heat 
energy of the atmosphere, we shall be able to find the energy supplied to the 
eddies, at least if the air be dry, by the aid of the theorem of Section 1 
given by Mr. W. H. Dines.* 

Now we can find the rate of change of intrinsic energy in terpis of the 



•X- t 



Q. J. Roy. Met, Soc.,' p. 188, July, 191^ 
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rate of eddying, by means of the equation for the diffusion of potential 
temperature, denoted by <l>, 



where p is pressure used as a measure of depth, and where |^ is a measure of 
turbulence, for which I have suggested the name " turbulivity." In quoting 
this equation from a previous paper it is necessary to make clear three points 
which I then^* left confused. For we are now examining small effects, and 
precision in details is essential to success. 

(i) When there are no eddies we are accustomed to compute the flow of 
heat or of water vapour across a plane from the flow of mass across the 
plane. As the effect of eddies is to be treated as additional, we should not 
include in it any flow due to the mean motion of mass across a plane, and 
accordingly D is here a differentiator which follows the mean motion of the 
fluid, in the sense that the effect of eddies is to cause equal masses to cross 
the co-ordinate surface from opposite sides. In the ease in which there is no 
convergence of air in the general circulation, D accordingly refers to a surface 
of fixed pressure, not to one of fixed height. Therefore, as ;p is the other 
independent, D/D^ may be replaced by curly djdt, indicating differentiation 
at constant pressure. 

(ii) Secondly, the equation would not hold if entropy per mass, <r, were 
substituted in place of (p, because entropy is increased by mixing. In fact, 
as for dry air, 

cl(T = ry^dlog^ (4-2) 

the equation for the diffusion of a takes the different form 

(iii) The eddy motion produces an atmosphere with an intricate structure ; 
but as the diversity is smoothed out by molecular diffusion .at constant 
pressure, the mean value of <I> for a secluded portion of air remains unchanged. 
Therefore in equation (4'1) we may regard <I> as a mean value over a portion 
of an isobaric surface. 

With this explanation let us return to (4*1), and deduce from it the rate 
of change of temperature. 

For dry air 

^ = 6>(2ij (4-4) 

where p» is the standard pressure adopted in defining ^. 

* *Roy. Soc. Proc.,' A, London, vol. 96, p. 10 (1919), 



Energy from and to Atmospheric Eddies, 361 

So, at constant pressure 

Insert this in (4*1), transfer the factor in p^'^^ to the other side of the 
equation, multiply by 7^ and integrate from the top to the bottom of the 
atmosphere, with respect to the element of mass per unit horizontal 
area clpjg 

Integrate the second member by parts 



g Jo 3^ 



*^ f: 



Pt^'^^ 9 Jo 3p 



The two terms on the right are both activities due to convection per unit 
horizontal area. They may be distinguished as the " boundary activity '' and 
the '' body activity " respectively. The boundary activity must vanish at the 
top of the atmosphere. The sign of the second integral in (4*7) determines 
whether, in the body of the atmosphere, heat is being transformed to 
eddying energy or vice versa, 

Now ^ is essentially positive and so is d(p^'^% since the integration is 
made downwards. Therefore (4*7) implies that heat becomes eddying energy 
if the potential temperature increases downwards. That is what we should 
expect. If there were no eddying, ^ would vanish and equation (4-7) would 
correctly indicate no activity. The atmosphere might then be unstable, but 
if so, it would be at rest in equilibrium — a hypothetical case merely. 

Consider the boundary activity in equation (4*7). Now, in the paper* 
referred to it is shown that ^/g^d^/dp is the amount of Z rising across 
surface per area per time, when Z is defined by the statement that <P is the 
amount of Z per unit mass of atmosphere, and the surface ought to have 
been specified as one crossed by equal masses from opposite sides. Thus 
equation (4*7) shows that the boundary convective activity is simply the rate 
at which eddies remove heat energy from the ground, vegetation, or sea per 
horizontal square centimetre. 

To correspond with the equation of continuity of energy in the form (llO) 
of Section I, we must integrate (4*7) over the large horizontal area A. When 
so integrated, the boundary activity per area becomes what we have formerly 
denoted by G-. x\nd the first member is 7^/7^, times the whole rate of increase 
in intrinsic energy, not arising from radiation, nor from dissipation of eddying 

^ *Eoy. Soc. Proc.,' A, London, vol. 96, p. 11 (1919) equation (16). 
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kinetic energy by molecular viscosity. So we obtain for the integral of the 
body activity, the following expression :- 

%p Sol 
7i; 3^ 



_a=-jr,ro-:^^j;i|#«- 



2» 



dK. 



(4-8) 



This expression can bo much simplified by using the following relations 



V 
dp 

da- 

\pl 



g^pG, where c is eddy conductivity 
hpd (gas equation), 

—gp (hydrostatic equation), 
^pcl log ^ = ^pcl log 6— hcl log^j, 

e 



# 



(definition of potential temp.). 



(4-9) 
(4-10) 

(411) 

(442) 
(413) 



The stages by which the transformation takes place may be as follows :- 






o'' ^ dp 



#. 



7^ i c-^ i -:^~-^,d^ 



dh L# ^^ 



3/^ I 



40— del log ^ 



'^ dor W dp __ g 



^ 9o- 



d PTT dli, (414) 



dh 7^ f jp] 

The last form is the simplest. When inserted in (4*8) the latter becomes 



yp dpi 



-G 



+ 



JL{\^1 



■ /-sj Cttv 

7^Ja oil 



WjcL. 



(415) 



In words this may be stated as follows :- — The average rate at which intrinsic 
and gravitional energy jointly are being transformed into the energy of eddies is 
equal per vohime to da/dhxcg/yp, that is to say to the product of the eddy con- 
ductivity into the accelercdion of gravity into the gradient of entro'py per mass 
iipward, divided by the thermal capacity per mass at constant pressure. 

Now there exist observations of c and of dcrfdh. So let us divide the 
activity per volume cg/jp x da/dh into the factors c and g/yp x da/dh and make 
a table of the latter under average conditions. For computing da/dh from 
observations, a convenient formula is 

dcr If 

dh'^d 



{7.1+^7} , 



(416) 



it applies accurately to dry air, and with quite small errors to clear moist air 
when yp has its value for dry air. So I have taken data for clear days where 
such data were available. It is not worth considering cloudy days until the 
theory can be arranged throughout for that purpose, and that would very 
greatly add to its complexity. 
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Table I. 
(Computed from Josef Reger s Analysis of the Lindenberg Observations,* 





clear days only.) 








Hour 
(0 = midniglifc). 


2. 


8. 


14. 


20. 


Layer, metres aboTe M.S.L. 
122 to 500 


gjyp d<ijcjh expressed in sec." 

10-4 X 
6 -52 3 -02 
1 -50 1 -95 
1 -83 1 -70 
2-51 1-20 


■2 equals 

the following 

-0-31 

0-76 

1-46 

1 -47 


2-54 
1 -07 
0-88 
0-20 


500 1000 

1000 1500 

1500 2000 



Table Il.f 
(England, S.E., days both clear and cloudy.) 



Height, 
km. 


gjjp 0(^1 vh in seo.'"'^ 


Height, 


gjyp vfrjvh in sec.~\ 


10-4 X 


km. 


10-4 X 


13-5 


4-42 


6 5 


1 -07 


12-5 


4-74 


5-5 


1-20 


11-5 


4-07 


4*5 


1-28^ 




10-5 


3-03 


3-5 


1 -47 


. •♦• 


9-5 


2-59 


2-5 


1-66 


r"+ 


8-5 


1-49 


1-5 


1-92^ 




7-5 


1-01 







J Doubtful, owing to formula for dry air having been used. 

These Tables show that on clear days in the troposphere gfyp2<T/dhiB of 
the order of 10"* sec.~^ on the average. 

Now the eddy conductivity at the height of a kilometre has been found by 
numerous observers to be of the order of 100 grm. sec."^ cm.~^. 

Taking these numbers it follows that Ci^/7^.9o-/3y^ == 10~^erg. sec."' cm~^. 
This is the convective activity per volume lost by the eddies. 

V. Joint affects of the Various Supplies. 

In the equation of continuity in the form (1*10) let us now insert the 
convective activity given by (4*15) together with Reynold's expression for 
the supply of energy from the mean wind (2*3). The result is 



dt 



MKSj 



+ /m(-^) YdhdA 



\dh/ J 



1 






^ 3cr " 

a oil 



dA. (5-1) 



■^ Vieweg und Sohn, Braunschweig, *Arbeiten des K. P. Aerowautischen Observa- 
toriums bei Lindenberg,' vol. 8, p. 247 (1912). 

t Computed from ' Characteristics of the Free Atmosphere,' by W. H. Dines, p. 62. 
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In reckoning the change of temperature due to convection, we took no 
account of the heat produced by the dissipation of eddy motion by molecular 
viscosity. So that dW /dt, as just given, ought to be decreased by the quantity 
(3'5) and a completer formulation is therefore : — 



t-^-' = m'-ijM'm ----Ef/i^* 



6vxjL» 



(5-2) 

How G. I. Taylor* put forward the hypothesis that the eddy viscosity fi = 
eddy conductivity c, and observation has so far tended to confirm this. Let 
us make this simplification. Then (5*2) becomes 

This is a very interesting equation. To take a special application of it, 
suppose that the large piece of atmosphere, throughout which the integral is 
taken, is initially at rest. Suppose next that an approaching depression sets 
it in motion. Will the winds be turbulent or not ? Suppose that they are 
very slightly turbulent, so that c and E' are just not zero, but so that 
2fjbc'rT^W l{pl^) is very small. The question is : will E' increase ? It will tend 
to increase if 

If the atmosphere has the ordinary value of g/jp'dcr/dh equal to say 
+ 10~*sec.~^ then the wind will become turbulent if the change of velocity 
with height exceeds ±y^IO~*sec.~^. That is to say, a change of about 
1 metre sec.""^ in 100 metres of height would be the critical value on an 
average clear day in the troposphere. Now pilot balloon ascentsf show such 
up-grades generally in the first 100 metres above ground and frequently up 
to a kilometre. But above a kilometre such up-grades are rather uncommon 
until the base of the stratosphere is reached, where again there is a rapid 
change-l But at this height, as Table II shows, g/yp'dcr/dh has increased 
considerably, so that a more rapid variation of wind, say 18 metres sec.""^ per 
kilometre of height, would be required to produce turbulence. The stability 
of the stratosphere is well shown by the large value of g/ypda/dh. We 
conclude then that the usual source of atmospheric eddies is in the first 
few hundred metres, and that they may sometimes be formed at the base of 
the stratosphere. 

The limit to E' will be set by the dissipation term 2fic'rr^B' / pP. 

^ *PliiL Trans.,' A, 1914. 

t * Manual of Meteorology,' by Sir Napier Shaw, Part IV, Ch. VI. 

J G. M. B. Dobson, ' Q. J. Met. Soc.,' Jan., 1920. 
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Stirring a dry atmosphere tends to bring the entropy-gradient, do-fdh, 
towards zero. So if the atmosphere is originally stable, equation (5'3) 
indicates that the more it is stirred the easier does it become to stir. 

The dying away of the surface wind at night has for some time been 
explained as an effect of radiation, which by increasing da/dh makes the lower 
layers too stable for eddies to exist in them. Equation (5*3) conforms to this 
explanation and gives to it an expression which is quantitative. Equation 
(5*3) also shows that, when the day wind is sufficiently strong (dv/dhy will 
overpower the term in dcr/dh, so that eddies will continue to exist and to 
supply momentum to the surface air, even on a clear night. That fits with 
observation. At Benson, for example, where the anemometer head is 
26 metres above the plain, the day wind usually dies away on a clear 
winter night if it is 5 metres/sec. but not if it is 11 metres/sec, provided 
that evenings on which the barometric gradient is changing rapidly are 
excluded. 

A theory like this one, which supposes the mean velocities to be in hori- 
zontal straight lines, can only fit with observation at a height above the 
, ground which is large compared with the irregularities of the surface. 




To make a proper test of equation (5-4) the following installation would 
be suitable : Two anemometers, one above the other, to determine dv/dh ; 
the lower should be at a height large compared with the local irregularities 
on the ground ; a difierential thermometer to read the difference of tem- 
perature between these two anemometer heads; and a pressure tube 
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anemometer placed at an intermediate level to record the range of the gusts, 
I am not aware of the existence of such an installation. 

Choosing an occasion at Benson^ when trees interfered as little as possible 
with the wind currents, I obtained the results shown in the figure. The 
temperatuie readings were taken with an aspirated thermometer. It is seen 

that the eqmroacJd of the ctcrves tor logio i— -^A and logio (?ct) eomdcles im 

"■ \jpO/i/ \6hl 

time at 16 h. 40 m. toith ci marked decrease in the range of the gusts 
shown hg the anemometer, A perfect agreement could not be expected 
because dvjdh was measured over the interval 1 m. to 26 m., while da/dh 
onlj over the interval 1 m. to 16 m. Again it would have been better if 
gustiness could have been measured at the average height, but trees pre- 
vented this. On this night the wind did not die awaj after sunset, but the 
current from 19 h. to 23 h. is probably the katabatic wind from the Chiltern 
Hills described by Mr. E. V. Newham in M. 0. Professional Notes No. 2. 



VI. the lrremTS%b%t%ty of Uonvection. 

We may think of this by the aid of the resulting increase of the total 
entropy S in an atmospheric column sfcandins upon a horizontal unit area. 
Since ^pfg is the mass in the short length of the column 



LuVj X. 



^d<T 



dp. 



dt g]o dt 
But (41) with (4*2) gives as an alternative form of (4 3) 

dt <lyV apl' 
Inserting this in (6*1) and integrating by parts 



(6-1) 



(6-2) 



^s^ 1 r» 1 a 

ell "' g Jo 



#qp\ op. 



d(r\ J ^ 
dp = 











9 



^ . da c3<i> 



^ dp ^%p 



-r-TT-- dp. (6 '3) 



The integrated term is simply the rate at which entropy is entering the 
column from the earth's surface. Eemembering - that da — ^^d\Q%^ the 
remaining integral may be written in alternative forms as follows: — 



cIS 



dt 



f 0(T 
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r« /aioo-#\2 

c — 5^ — dk 
G \ oh 



1 f^ 'd(T\^ 

c (-^-7 dh, -.6*4) 
1 \oh 



jp 



Q 



This quantity is essentially positive as we should expect, since the 
entropy cannot decrease. It shows that the rate of increase of entropy 
depends upon the square of da/dh, whereas the rate at which eddying 
energy is converted iuto heat depends upon the first power of the same, 
see (4*15). Thus the nearer the gradient of entropy approaches to zero the 
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more nearly will the eddies approach to^ perfectly reversible engines. If 
the eddy conductivity c were confined to a range of height in which da/dh 
were constant, so that^ we might remove the integral signs, then dividing the 
last form of (6*4) by the last of (4'14) it wonld follow that 

increase of entropy __ 1 3(t 



energy supplied from eddies g dh' 
Multiplying by the temperature of the layer we should find 



(6-5) 



energy made unavailable __0 (^(r .«./>a 

energy supplied from eddies g dh' 

This fraction is a pure number, which has an average value of about 
one-third in the troposphere. We see from (4-16) that in an isothermal 
atmosphere it would have the value unity. In an adiabatic atmosphere it 
is zero. 

VII. Eddy Thermodynamics. (Kevised May 11, 1920.) 
The direct eddy stresses are, according to 0. Keynolds 



XX ~ -^pv^'v^ ' \ yy =z ^pv^ ' v^ ' ; Jih = —pv^'v^\ (71) 

when tractions are reckoned positive, and bars denote mean values, and 
dashes denote deviations from the mean. 

Eollowing custom in the theory of viscous fluids, we could define an " eddy 
pressure," ^, as the negative mean of the direct tractions 

<^Q := — i(OT + 2/^ + M). (7*2) 

From (7-1) and (7*2) it follows that 

* 

^ = t . H^/^ + %'^ + ^h5 (7'3) 

P 

This resembles the " gas equation " '^^ = b0 and suggests that we might 
tentatively introduce the idea of *' eddy-heat per mass," defined thus 



" eddy-heat per mass " = | ('^x '^ + % '^ -f Vj^ '^) = B say, (7-4) 

so that p=:^p®. (7-5) 

Next taking a suggestion from the adiabatic cooling of gases, let its 
enquire whether the eddy -hectt-jper -mass would diminish if the air Gontairting 
the eddies were suddenly expanded. In this connection we must note that 
the energy per mass O, with which we are concerned, is not acoustical, 
and so moves with, not through, the fluid (Helmholtz, Kelvin). It is true 
that @ as defined by (7*4) would include some acoustical energy, but let us 

2 D 2 
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suppose that sufficient time has elapsed to allow this to disperse itself wifch 
the speed of sound. 

To obtain some light on the question, let us imagine a region containing 
a lively eddy motion accompanied by a mean-expansion. Eeynolds' investi- 
gation relates to an incompressible fluid, so we must make a fresh start. 

Consider a plane moving so that on the average as much mass crosses it 
from one side as from the other. For simplicity let the plane be set at 
right angles to v^. The molar kinetic energy in the fluid at any point is 
2P('yx^4-'yY^ + 'yH^) per volume. And not being acoustical, it is moving with 
the velocity of the fluid %• As the plane is moving at the speed %, the 
volume of fluid crossing an area dA. of the plane in unit time is (% •— '^x ) '^A. 
So the mean rate at which molar kinetic energy crosses the plane is 

i|ip(^xH^YH^H')(^x-'^x)^A. (7-6) 

Now introduce dashes for variations from the mean so that 

v^^ v^-\- v^^ ] p :=z p + p I and the like. (7*7) 

Then (7'6) becomes, on squaring out, 

+ t;/2 + ^Y'^ + ^H>x'^A. (7-8) 

Now the mean values of the rapidly varying quantities, which are indicated 
by dashes, will approximately vanish when taken over a sufficiently large 
area. So also will vanish the mean of any dashed into any barred quantity, 
since the barred quantities vary slowly from point to point. Thus we may 
pick out from (7*8) the following terms which vanish : — 

\lhp{v^' + v,''+Vn'WdA = 0. (7-9) 

Again, since an increase of density is as likely to be accompanied by 
a positive as by a negative %', so the following terms in pv^' vanish : — 

~||/ {v^'-\~v,' + 'V^'WdA == 0, (7-10) 

which is equivalent to assuming the absence of convection of mass. 
In order to divide our difficulties, let us at this stage assume that 

Tlip(%'' + %'M-^H'^)%'^A = 0. (711) 

This assumption implies that the eddying kinetic energy per volume is not 
correllated with v^\ so that portions of high energy are not moving specially 
to either side of the plane. That is to say, we have assumed no diffusion of 
eddying energy, from where it is intense to where it is feeble. We shall 
return to the question of diffusion in the next section. 
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Let (7*9), (7'10), (7*11) be assumed, then (7*8) simplifies to 

Yli^xP'^x^'^x' + Vy p'i^Y't^x' ^'^j,pv^' v^') dA, (7-12) 

Now, in conformity with Eeynolds' statement of the eddy stresses, when 
the bar replaces ~ |( )dA as an averaging operator, put to complete (7*1). 



xy =: —pVx^'Vy'; yh =^ —pVy'vi/] te == —pVjj^'v^', (7*13) 

Then the mean rate at which molar kinetic energy is crossing unit area of a 
plane set at right angles to v^ and moving with the fluid at its centre, 
becomes from (7*12) 

— '^x • XX— Vy . xy—v^, xli, (7*14) 

Now, if we added to the plane, which we may imagine as square, five other 
planes, also moving with the fluid to form at one instant a unit cube ; then 
we see that the rate of increase of molar kinetic energy in the moving, 
distorting, swelling cube will instantaneously be 

■^ {v^ . XX ^- Vy . xy + '^H . xh) 4- ^ ('^x • yx + Vy yy + "o^ yh) 

n --^ ^-^ '-^ 
H- ^ {v^ . hx 4- Vy % + '^H ^^^^)' (7*15) 

Since the mass of the cube is p and does not change as the cube swells, the 
expression (7*16) must be equal to 

P^^ihiv.'+v.J + v^^)}, (7-16) 

where I) denotes a differentiation following the mean motion. This gives the 
decrease in the total molar kinetic energy. We want the part of it associated 
. with the eddies, separated from that associated with the mean wind. But 
the changes which we have been considering are between the cube and fluid 
touching it, so that they concern only the relative velocity of the cube and 
its surroundings. Therefore without loss of generality, we may take v^, Vy, v^ 
as vanishing at the centre of the cube. Then the energy associated with the 
mean motion vanishes, leaving only the part in ©. So (7-15) when equated 
to (716) yields 



XX r\ ~J~ Xcf r\ ~I~ X/v r\ 

OX ^ ox ax 

-^ dvy^ ^^ dvY ^ B'^H 



D 



hx -T— + hy -7— + hh -?~ 
oh "^ oh oh 



y 
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The first member is Eeyiiolds' Activity, formed of the products of the eddy- 
stresses into the corresponding components of mean strain. We have arrived 
at it by a proof which differs from Reynolds' in that we liave : — 

(i) Supposed the fluid to be compressible. 

(ii) Done away with the rigid envelope at which he made all the velocities 
to vanish, and have replaced it by certain assumptions as to the 
approximate vanishing of various means. Just as ordinary thermo- 
dynamics is correct because the number of molecules is enormous, 
so eddy-thermodynamics is an approximation which improves as 
the number of eddies increases. 

In so far as Eeynolds' Activity depends upon distortion, we have already 
dealt with it in Sections II and V above. Let us now suppose that the 
rateS"Of-mean-shearing are zero, and let us examine the remaining part due 
to the expansion. Accordingly, (7*17) will be limited here to 

^dv^ ^dVy f^9^H T)B /hrlQX 

^x ^'^ 3z/ m ^ ^^^ ^ ^ 



Now observations by G. I. Taylor and by the author* have shown that, 

except in the first few metres above the surface, xx = yy ^ hh, at least 
roughly. In other words, there is a tendency to equi-partiti on of energy 
between the three components of velocity. For simplicity, let us assume 
this to be the case, so that from (7*2) 

XX = yy = hlh = — ^. (7'19) 

Then from (718) 

(7-20) 



l^dx dy dh 




1)0 

-p m 


But by the continuity of mass 






dv^ dv^r dv^ 
dx ' dy ' dh 


Dt 


(log p), 


therefore from (7-20) and (7-21) 






S^j3,Iog/> 







(7'21) 



(7-22) 



03?, using (7 '5) and transposing 



Dt 



{|log/9~~ioge}==:0, (7-23) 



or ^~~. is constant following the mean motion. 

That is exactly like the adiabatic equation for a gas having a tem- 

^- ' Phil. Trans.; A, 1920. 
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perature @ and a ratio of specific heats equal to 5/3 — like mercury vapour, 
krypton, and other monatomic gases. 

One is tempted to call % the " eddy-temperature," but it is perhaps best to 
avoid this term, for @ has different dimensions, and so cannot be expressed in 
degrees on a thermometer. It is correctly called '' eddy-heat-per-mass/' 

We are thus led to the conception of a '* potential eddy-heat-per-mass," 

defined as 

@ {pilpfl^ (7-25) 

where pi is some standard density. 

In arriving at (7'14), we ignored the fact that the molecules do not all 
move with the velocity %, %, Vh, which we attribute to the fluid at a point. 
But these molecular motions are the cause of the ordinary pressure p. By a 
method, closely analogous to that by which (7*14) and (7*15) are obtained, it 
can be shown that (7'14) should be increased by —v^;p and (7'15) by 

r\ P) r\ 

^(^'xi^) + ^ ('^yi^) + ^ ('^Hi?) + terms depending on molecular viscosity. 

Odb (jlj (Jib 

These in turn lead to the lowering of the ordinary temperature at a point. 
But all that is well understood, and our simplest plan is to leave out f from 
all the equations, and to expect no result relating to ordinary cooling. 

Similarly, with regard to gravity. That has been left out. If it had been 
included, we should have had to consider the mean flux of gravitational 
energy over a plane moving with the mean motion. The flux will vanish if 



p'^h' - 0, (7-26) 

that is to say, if convection can be neglected. But having already dealt 
with convection in Sections IV and V, we may omit it here. Thus the 
increase of molar-kinetic and of gravitational energy in the cube is still given 
by (7*15). But, in place of (7*16), we shall have 

P^{t + H^x^ + ^r^ + ^H^)}, (7-27) 

where i/r is the gravity potential. 

It is not now permissible to suppose v^ zero at the centre of the moving 
cube, since v^ must be relative to the attracting earth. So the first member 
of (7*17) must be increased by 

while the second member of the same equation is increased by 

P§i{i^ + V"n'). (7-29) 
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But these two increases balance one another, on account of the dynamical 
equation. Therefore (7*17) is correct also when gravity is present. 

The foregoing theory, which states that the eddy-heat-per-mass is 
diminished by a mean expansion, may be expected to have some application 
to the following natural phenomena : (i) turbulent winds flowing up 
mountain sides ; (ii) numerous small eddies carried together upwards in some 
larger circulation, such as that made visible by cumulus or thunder cloud. 

But on account of the assumption (7*11) the theory does not apply to 
the case in which the small eddies wander individually from crowded to 
roomy regions. 

VIII. Diffusion of Eddying Energy, (Ee vised May 11, 1920.) 

Lieut.-Col. A. Ogilvie* describes flying over the Mle in a wind of 13 metres/ 
sec. : 'Howards 10 o'clock or so the river appeared like a sheet of glass, but 
a little later, when the sun began to heat the ground, one saw ripples coursing 
over the surface of the river," and after half-an-hour or so the gusts reached 
the aviators at a height of about 2/3 km. 

Capt. 0. K. M. Douglasf also writes of '' rising turbulence " as of some- 
thing commonly experienced. These observations fit in with the theory of 
Section V to the effect that turbulence usually originates near the ground. 
According to What law then does it diffuse ? Two different processes may go 
on at the same time. Eddies may wander individually from places where 
they are crowded to places where they are rare. We may speak of this as a 
true diffusion of eddy-heat-per-mass. Or large whirls, such as those indicated 
by cumuli, may lift several' cubic kilometres of air containing thousands 
of smaller eddies. This might be called a large-eddy-diSusion of small-eddy- 
heat-per-mass. Let us discuss this first. 

Smallsccde' Eddy 'heat Carried up hj Larger Eddies. 
Can it, for instance, be fitted into the following form of diffusion equation ? 

where ^i is some unknown quantity which is not zero when d^/c^p = 0. For 
this to be a possible form ^j^ must satisfy the following conditions.^ 

1. ;j^ must be carried along with the moving air. (8'2) 

2. X i^^st be unaltered by the expansion or contraction of the 

air with which it moves. (8-3) 

^ ' Aeronautical Journal,' p. 415, July, 1919. 

t ^ J. Scott. Met. Soc.,' vol. 17. 

I Taken with revision from * Roy. Soc. Proc.,' A, vol. 96, pp. 10, 11 (1919). 
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And if dm be an element of mass and Ix^lm be taken over a definite portion 
of air, then 

3. J%<im must be unchanged by internal rearrangement of the 

portion, at constant pressure. (8'4) 

4 Ix^'^'^ must be unchanged by delay. (8*5) 

May we put for % the " potential-eddy-heat-per-mass " discussed in 
Section Til ? That would satisfy (8*2), (8'3X (8-4). Condition (8*5) is not 
at all satisfied unless we treat the diffusion-process as supplemental to the 
changes indicated in (6'2) but with this understanding (8'5), and therefore 
also (8'1), would be satisfied if 

^ = potential-eddy-heat-per-mass = B ( ^ | , 

where % is the '* eddy-heat-per-mass " and pi is the standard density adopted 
for the purposes of definition. 

The true diffusion of eddy-heat bears the same kind of relation to the 
convection of small eddies with cumuli as molecular conduction of ordinary 
heat bears to eddy-conduction of the same. We have already encountered in 
(7'11) an expression giving the mean rate at which energy diffuses across 
unit area of a plane perpendicular to the ^— axis, in the form 

but it is required to transform this into another applicable to observations. 
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